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On fractals, spectral functions such as heat kernels and zeta functions exhibit novel features, very 
different from their behaviour on regular smooth manifolds, and these can have important physical 
consequences for both classical and quantum physics in systems having fractal properties. 

I. INTRODUCTION 

Stuart Dowker has made many fundamental contributions to the modern understanding of heat kernels and zeta 
functions, and their role in theoretical physics, ranging from quantum field theory on curved space such as De Sitter 
and anti De Sitter space [1 , to vacuum energy and Casimir physics [2] , to finite temperature quantum field theory [3] . 
A common theme in his work is the great utility of these spectral functions in encoding physical properties of quantum 
systems confined to interesting curved manifolds, such as spheres, cones and orbifolds, often highlighting surprising 
and deeply beautiful connections to number theory. Having personally learned so much from his papers, it is a great 
pleasure to contribute this paper in honour of his 75th birthday, and to point out some novel features that occur when 
such smooth manifolds are generalized to fractals. The most direct physical application is to classical and quantum 
transport in fractal or aperiodic media, which exhibit many novel phenomena fHUO) . A more ambitious application is 
to quantum gravity, both for two-dimensional models and in four dimensions. Recently, it has been shown that 
a wide variety of different approaches to quantum gravity, such as asymptotic safety and the renormalization group, 
causal dynamical triangulation and Euclidean dynamical triangulation, all show hints of some fractal-like features of 
quantum gravity at short distance scales [T2] EH H5HT7] . 

Heat kernels and zeta functions are well known in terms of diffusion processes and path integrals. This is very 
well understood in flat space, but becomes more interesting and non-trivial when particles and fields reside on curved 
spaces or couple to gauge fields. The key philosophical idea is to relate spectral information of the propagating 
fields or particles to the geometric properties of the manifold [or fractal]. It appears that the zeta function is a 
particularly convenient and efficient way to encode this spectral information, while the heat kernel enters the physical 
story primarily because of its direct relation to the quantum path integral, and to the statistical mechanical partition 
function. In a very real sense, we use the propagating fields and particles to probe the spatial geometry. An early 
example of this, to be elaborated further below, is Lorentz's question: why does the Jeans law only depend on the 
volume of the region in which a gas, at thermal equilibrium, is contained? This provided physical motivation for 
Mark Kac's famous question: "Can one hear the shape of a drum?" . Perhaps the most primitive and basic property 
of a manifold is its dimension. As is familiar from physics, the behaviour of a system, especially a quantum system, 
is profoundly dependent on its dimension. Fractals, having fractional and anomalous dimensions, are a natural place 
to look for exotic physics and mathematics. 

The mathematical study of diffusion on fractals was motivated by results from the physics community, where the 
concepts of spectral dimension and walk dimension were developed [TBH2TJ] . In the late 1980s a general mathematical 
theory was developed for diffusion and harmonic analysis on fractals [2TH2rj] . More recently, results have been proved 
concerning existence and uniqueness of Brownian motion and the Laplacian on a wide variety of fractal spaces. There 
have also been recent advances proving existence and meromorphic continuations of the zeta function on certain 
fractals, with rigorous results and explicit expressions for the Sierpinski gasket [27, 28J. A general picture is emerging 
in which the zeta function for highly symmetric fractals has complex poles, and these result in log periodic oscillations 
in the small time asymptotics of the associated heat kernel trace. Log-periodic oscillations have a long history in 
physics: in the theory of phase transitions [25], the renormalization group [30], Levy flights and fractals .31], and 
generally in systems with a discrete scaling property [32 . In this paper, dedicated to Stuart Dowker, I first review 
the familiar properties of zeta functions and heat kernels on regular smooth manifolds, and then I show how some of 
these change when generalized to fractals. 



* Invited contribution to the JPhysA Special Issue in honour of J. S. Dowker's 75th birthday. 
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II. ZETA FUNCTION AND HEAT KERNEL: BASIC PROPERTIES ON SMOOTH MANIFOLDS 



Before coming to spectral functions on fractals, I review the most relevant basic properties of the zeta function 
and heat kernel on smooth manifolds p5 38J . The later work on fractals will refer to the spectral properties of the 
Laplacian operator, but we can introduce the main spectral functions for a more general self-adjoint operator H . The 
zeta function is defined as 

= tr^7 = E^ • W 

The sum over the spectrum {A} converges for the real part of s large enough, but we are particularly interested in 
the question of whether or not £(s) has a meromorphic continuation throughout some further region of the complex 
s plane. For example, in Casimir physics, where the relevant operator is a wave operator, with eigenvalues w 2 , we 
are interested in making physical sense of the formally divergent sum for the zero point energy: Yl \^ lL0 ~ C( — 1/2). 
In computing the quantum field theoretic effective action, we need the logarithm of the determinant of the operator 
H [perhaps a Schrodinger operator, or a Klein-Gordon operator, or a Dirac operator, . . .], which is formally related 
to the derivative of the zeta function at s = 0, lndeti? ~ — £'(0), as was first defined by Dowker and Critchley pQ. 
Making sense of such formal expressions is the physical challenge, and the zeta function plays a privileged role. 

Another important spectral function, which enters because of its relation to the Euclidean path integral and the 
thermodynamic partition function, is the heat kernel trace: 

K(t) =tre- Ht = ^e- At . (2) 

A 

The heat kernel trace and the zeta function are transforms of one another: 

1 f°° dt 

C(s) = — / -t s K(t) , (3) 

r(s) Jo t 



i 
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K(t) = — I dst- s T(s)as) ■ (4) 



where the contour for the inverse Mellin transform in Q is determined by the analytic properties of the zeta function. 



A. A simple illustrative example: the line segment 

A simple example, which still captures much of the essence of the situation on smooth manifolds, is the Laplacian 
on the interval of length ir. The eigenvalues of the Laplacian are X n — n 2 , with n a non-negative integer starting at 
n = for Neumann boundary conditions, but starting at n = 1 for Dirichlet boundary conditions. Then the heat 
kernel trace is [— for Dirchlet, and + for Neumann]: 



K(t) = \ jr e-"*=F 



1 



2 '2 

Tt — — OO 

I / 00 1 



2V i ^ 2 

n— — oo 



e-^T 

(5) 

The factor n in the numerator of the first term is the volume of the one-dimensional interval [cf ^ below] . The second 
term is a boundary contribution, whose sign depends on the type of boundary condition. The zeta function is also 
very simple. For Dirichlet boundary conditions (or, separating the zero mode, a subject to which Stuart Dowker has 
made important contributions [39] ) . we see that the zeta function is expressed in terms of the well-known Riemann 
zeta function: 

COO = E 7^1 = (n(2s) (6) 

n=l V ' 
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This zeta function has a well-known analytic continuation into the complex s plane, being meromorphic with a single 
simple pole at s = 1/2. From the transform expression Q we see that the contour must be taken to the right of 
the line Re(s) — 1/2, and when deformed around the poles of £(s) and T(s) we pick up the two terms of the small 
t expansion of K(t) in ^ [note that the simple poles of T(s) at the negative integers do not contribute because 
(r(— 2n) = 0]. Thus, as can be understood from the transforms plfl]), we learn an important lesson: 



poles of £(s) < — > small t asymptotics of K(t) (7) 

B. Smooth Riemannian manifolds: the Weyl expansion 

It is straightforward to generalize the previous example to a <i-dimensional hypercube, and from there to any smooth 
Riemannian manifold. This leads to the fundamental result on smooth manifolds, known as the Weyl expansion 
[351 00], that on a d-dimensional manifold the heat kernel trace has asymptotic small t behaviour: 

K(t)~ V , 2 -a S _ t + ... (8) 
{4irt) d / 2 (47rf) — 

where V is the volume, S is the boundary volume, and a is a numerical constant determined by the boundary 
conditions. The dots refer to higher order terms depending on further geometric properties of the manifold such 
as curvature. The leading term of the expansion ^ relates spectral properties of the Laplacian to basic geometric 
properties of the manifold: its dimension d, and its volume V. We can understand the small t asymptotics in (8| 
from Weyl's result [41] concerning the large n growth of the eigenvalues, A„ ~ n 2 / d , for a self-adjoint second-order 
differential operator on a (i-dimensional smooth manifold: 

K(t) J dne-*"' (9) 

71 

Correspondingly, 

<M~E(^~0>(f) do) 



which has a simple pole at s = d/2, which, using the transforms in (J3j [4J) , is directly related to the leading small t 
asymptotics in 

Two natural questions to ask are: on a fractal, is there such a Weyl expansion, and if so, what plays the role of the 
dimension d, and the volume VI These questions are addressed below in Sections 4 and 5. 

C. Heat-kernel propagator and off-diagonal Green's function 

Another physically interesting spectral function is the off-diagonal heat kernel, 

K{x,v-t) = ( V \e~ Ht \x) = P t {x,v) , (11) 

which is the propagator P t (x,y) for the associated heat equation 

d t u = -Hu (12) 

In probabilistic terms, Pt(x,y) is the probability for Brownian motion or diffusion to propagate from x to y in time 
t. In d-dimensional flat Euclidean space, the familiar result is 



K ( X >V> t )= eX P 



\x - y{" 



4t 



(13) 



which is clearly consistent with ([8[). For a more general manifold, this generalizes [35j[42]: on a geodesically complete 
Riemannian manifold with positive Ricci curvature the small t asymptotic behaviour of this propagator is 



K(x,y;t) - 1 r exp 

V{X, y/t) 



c P 2 (x,y) 



(14) 



4 



where p(x, y) is the geodesic distance between x and y, c is some positive real constant, V(x, y/t) is the volume of a 
ball of radius yfi centered at x, and we assume that t is small enough that there is a unique geodesic connecting x 
and y. 

A closely related quantity is the off-diagonal Green's function of the Laplacian: 

i r°° 

G(x,y) = (y\-\x) = J dtK(x,y;t) . (15) 
On flat space we have the familiar expression 

G(x,y) = < 4 v d/2 \ x - y \d-2 , a T L ^ 
[±]n\x-y\ , d = 2 

Similarly, on a smooth Riemannian manifold we have the short distance behaviour 

G(x,y)~{wyW^ > f^ 2 (17) 
[\np(x,y) , d = 2 

These results highlight the importance in physics of the critical dimension d = 2, leading to many interesting physical 
phenomena in planar systems. For propagation on a fractal, we will see below how this critical dimension generalizes 
in an interesting way. 

D. Other spectral functions: Schrodinger kernel trace, Poisson kernel trace 

To study both classical and quantum transport on a manifold we often require information about the Schrodinger 
kernel trace 

S(t) =tve- lHt = ^e- lXt , (18) 

A 

which is clearly related to the heat kernel trace K(t), but has completely different convergence properties, being highly 
oscillatory. From the resolvent Green's functions 

G ± (E)^ W -^— , e ^0+ , (19) 



we find the density of states 



Further, defining the S-matrix as 



a{E) = tv 5{H - E) = -lmG+(E) (20) 

7T 



we can express the density of states as 

<r(E) = TT^^tr In S = J-^- In det S (22) 
v ' 2nidE 2m dE v ; 

On the other hand, the net transmission probability [for one-dimensional transmission] is given by 

\T(E)f = l/dct(G + (E)G-(E)) (23) 

Thus, both the density of states and the transmission probability are expressed in terms of the determinant of the 
resolvent Green's functions, which have the "proper-time" form in terms of the Schrodinger kernel trace: 

r°° a+ 

In det G±(f?)=trInG± (£?) = - / — e^^* tr e* iHt (24) 
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Finally, it is clear that other spectral functions, such as the Poisson kernel trace 

P(i)=tre-^ t = ^]e- VX * (25) 

A 

can be expressed as a transform of the heat kernel trace 




Similarly, the Poisson kernel trace V{t) can be expressed in terms of the zeta function, and so its small t asymptotics 
is also governed by the pole structure of the zeta function. 



III. SOME FRACTAL BASICS: HAUSDORFF DIMENSION, SPECTRAL DIMENSION AND WALK 

DIMENSION 

On fractals, the behaviour of these spectral functions is very different. The zeta function £(s) develops complex 
poles in the complex s plane, and these lead to log periodic oscillatory behaviour in the small t asymptotics of the 
heat kernel trace. In the next section we describe these results using an explicit example of the class of diamond 
fractals fi3l 03] , but we note that these features are quite generic to a wide class of fractal systems 125] . Before 
discussing this, we recall various natural dimensions associated with fractals. 

1. Hausdorff-Besicovitch dimension: this dimension, dh, usually referred to simply as the Hausdorff [or fractal] 
dimension, is the most familiar dimension associated with a fractal, and refers to the spatial scaling properties 
of the fractal: 

In V(r) 

d h = lim ^2 , (27) 
r->0 in r 

where V(r) is the volume of the fractal at length scale r. For a smooth manifold dh reduces to the usual 
integer-valued dimension d, but for a fractal dh is generically non-integer. For example, for the Sierpinski gasket 
shown in Figure [I] the Hausdorff dimension is dh — In 3/ In 2 ps 1.585, because on each iteration we halve the 
length of each interval but the number of intervals grows by a factor of 3. 




FIG. 1: Sierpinksi gasket embedded in two dimensions. This fractal has Hausdorff dimension dh = In 3/ In 2 w 1.585, and 
spectral dimension d 3 = 2 In3/ln5 w 1.365. 

2. Spectral dimension: this dimension, d s , refers to the scaling properties of the eigenvalues of the Laplacian 
defined on the fractal [I5H2D] . It is therefore tied directly to the propagation of heat via the heat equation 
dfU = Ait. A simple [but not quite accurate, as we discuss in the next section] way to introduce it is through 
the small t asymptotics of the heat kernel trace: 

d In K(t) 

d s = -2 lim — - w . (28) 
t^o dint y J 

For a smooth manifold, d s coincides with the integer- valued dimension d, and also with the Hausdorff dimension: 
d s = dh = d. However, on a fractal, d s is generically different from dh, and also is generically non-integer valued. 
For example, on the Sierpinski gasket shown in Figure [I] d s — 2 In 3/ In 5 « 1.365. 
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3. Walk dimension: this dimension, d w , is not really an independent dimension, but is a ratio of dimensions, and 
has the physical meaning of a diffusion index |46) . It encodes the anomalous diffusion that occurs on a fractal, 
where the usual Einstein relation for Brownian motion is generalized as [TSl EI] : 

(r 2 (t))~t — > (r 2 {t)) ~* 2 /<*» . (29) 

Generally, d w > 2 on a fractal, so that diffusion spreads more slowly on a fractal than on a smooth manifold; 
it is called "sub-diffusive" or "sub-Gaussian". In fact, d w is not independent of the Hausdorff and spectral 
dimensions. The anomalous diffusion implies that the probability to diffuse a distance r in time t is a scaling 
function: p(r,t) = f(r dl ° /t). But the normalization of total probability involves integrating over the fractal, 
which necessarily introduces the Hausdorff dimension, written symbolically as: J d dh rp(r,t) = 1. This then 
implies the scaling behaviour 

Comparing with the definition of the spectral dimension we see that 

d w = ^ . (31) 

For example, both numerical and analytical work shows that diffusion on the Sierpinski gasket occurs with 
d w = In 5/ In 2 s» 2.32, in agreement with the values of dh and d s mentioned above. 



IV. COMPLEX POLES OF THE ZETA FUNCTION AND LOG PERIODIC OSCILLATIONS OF THE 

HEAT KERNEL TRACE 

In fact, for a symmetric fractal, the small t behavior of the heat kernel trace is more involved than the leading 
term K(t) ~ l/t d / 2 in for a regular smooth manifold. Instead, both the spectral dimension d s and the Hausdorff 
dimension dh play a role in the leading small t behavior. For very general fractals, the heat kernel trace K(t) is 
bounded above and below by such a functional form: there are constants c\ and C2 such that [231 124) 

^ < m < ^ . (32) 

Mathematicians write this in short-hand as 

K{t) x ^ . (33) 

For many fractals, those with exact self-similarity [including for example the Sierpinski and diamond fractals], the 
result is in fact far more interesting than this, as K(t) exhibits characteristic log periodic oscillations between these 
upper and lower limits: 



, 2vr , 

1 + a cos - — - — - In t 



(34) 



K d w \nl 

for some real constants c, a and 0. Recall that the walk dimension d w depends on the ratio of the Hausdorff and 



spectral dimensions, as in (31). 

This type of log periodic oscillatory behavior is show in Figure [2j Note that the log periodic oscillation depends on 
the walk dimension d w and on a spatial-decimation factor I, which for the Sierpinski gasket is I = 2 [because in each 
iteration of the fractal each unit is divided by 2] . Rewriting ([34]) as 



ca I e~* v 



K ®~#Q2 + Tl ™ + «- ™ )+ -"- • (35) 



I 2 -l~d m lni f 2 



and recalling the transforms (|3j [4J> that relate K(t) to the zeta function £(s), this form of K(t) can be identified with 
a tower of complex poles of the corresponding zeta function: 

d s 2nmi 

Sm = Y + dZhd ' mEZ (36) 
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K(t)//q ea ding(t) 

1.0 




0.25 



FIG. 2: The log periodic oscill ations, at small t, for the heat kernel trace K (t) on a fractal, normalized relative to the leading 
term -Ki ca din g (i) = c/t da ^ 2 in (34 351. The solid [blue] curve is exact; the dashed [red] curve is the first two terms in the 



approximate expression ( 35 I 



as shown in Figure [3j It turns out that the contribution of the higher complex poles in the tower is exponentially 
suppressed by the Gamma function factors coming from Q, so that (35) gives the leading terms and in fact the 



second, log periodic, term is much smaller in magnitude than the first term, as can be seen clearly in Figure [2] 

Thus, d s describes the averaged leading behaviour of K(t) in the small t limit, while d w is needed to characterize 
the log periodic oscillations. Another definition of the spectral dimension is as the [largest] real part of the complex 
poles of the zeta function: 

d s = 2maxRe(pole of C( s )) ■ (37) 

We qualify this definition by saying "max" because there may be other complex poles, but with smaller real part. 
Furthermore, we see that the identification ([7| of the poles of the zeta function with the small t asymptotics of the 
heat kernel trace generalizes in an interesting way to: 

complex poles of £(s) < — > log periodic oscillations of K(t) (38) 



A. Illustrative example: diamond fractals 

These properties of spectral functions on a fractal can be illustrated in explicit detail using the class of fractals 
known as "diamond fractals", as illustrated in Figures [4] and [5j The diamond fractals are defined as iterations of 
graphs [43l |44j. At each step n of the iteration, we characterize the diamond fractal by its total length L n , the 
number of sites N n , and a physical diffusion time T n . Scaling of these dimensionless quantities allows to define the 
corresponding Hausdorff dh, spectral d s , and walk d w dimensions according to 

A lnNn A lnNn A lnT « (1Q\ 

dh = ^r n ' ds=2 wf; ' dw = ^L~ n ' (39) 

where the limit n — > oo is understood. These three dimensions are thus related by d s — 2dh/d w as in ( |31[ ). The 
diamond fractals are characterized by two integers, one of which describes how many times each interval is branched, 
and the other of which describes the decimation of each interval. Figure [4] shows the first two iterations of the diamonds 




Ll 

d w In/ 

dj2 



FIG. 3: Sketch of the complex pole structure of the zeta function for a fractal. Note the appearance of a tower of complex 
poles, whose real part defines the spectral dimension d s /2. On certain fractals there may also be a tower of complex poles 
with Re(s) = [47] . and there is a meromorphic extension of for all Re(s) > — e. Further results on the meromorphic 
continuation have been obtained by Kajino |48) . 



^4,2) Dq 2 and Dq 3. The further iteration of D4 2 is shown in Figure |5| w hich makes more clear the origin of the 
name "diamond" fractal. The associated dimensions are listed in Table IV A Notice that by choosing the decimation 
and branching factors appropriately the spectral dimension can be tuned as desired. Importantly, d s can be made 
either smaller than or larger than [or equal to] the critical value 2. Also notice that the diamond fractals are special 
because they all have walk dimension d w = 2, which means that d s = d h : the spectral and Hausdorff dimensions 
coincide. This is a special feature, but it does not affect the relevant properties of the spectral zeta function or heat 
kernel trace. 





dh 




d s = 2dh/d w 


I 


£>4,2 


2 


2 


2 


2 


A>, 2 


In 6/ In 2 


2 


In 6/ In 2 


2 


£>6,3 


In 6/ In 3 


2 


In 6/ In 3 


3 


Sierpinski 


In3/ln2 


In 5/ In 2 


21n3/ln5 


2 



TABLE I: Fractal dimensions and size scaling factor I for diamond fractals, compared with those for the Sierpinski gasket. For 
Dq ; 2, the spectral dimension is d s ~ 2.58, and for Dg^, d s ~ 1.63. 



The heat kernel trace for the Laplacian on a diamond fractal can be obtained by noticing that the spectrum is 
the union of two sets of eigenvalues [121 |33]- One set is composed of the non degenerate eigenvalues 7r 2 fc 2 , (for 
k = 1, 2, . . .). This corresponds to the spectrum of the diffusion equation denned on a finite one-dimensional interval 
of unit length, with Dirichlet boundary conditions. The second ensemble contains iterated eigenvalues, / K 2 k 2 Lt w , 
obtained by rescaling dimensionless length L n and time T n at each iteration n according to Lf™ = T n , given in ( 39 1 . 
To proceed further, we use the explicit scaling of the length L n = l n upon iteration. These iterated eigenvalues have 
an exponentially large degeneracy given, at each step, by BLff = B (l dh ^ , where B = (7 dh ~ 1 — 1) is known as the 
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FIG. 5: Further iteration of the diamond fractal 1)4,2- 



branching factor of the fractal, and the integer l dh is the number of links into which a given link is divided. The 
exponential growth of the degeneracy plays a crucial role in our analysis. Then the diamond heat kernel trace Ko(t) 
is the sum of contributions of the two sets of eigenvalues: 



00 00 00 

(,2 2 



K D {t) = J2 e + B J2 indh E e_fe • ( 4 °) 



k=l n=0 fe=l 

Correspondingly, the associated zeta function, Cd( s )j is 



n=0 



where (j{(2s) is the Riemann zeta function. Note that the zeta function has a factor involving a geometric series; a 
very similar structure arises for the zeta function of the Sierpinski gasket 27, 28]. The geometric series factor in (41 ) 
implies that the diamond zeta function Cd( s ) nas a tower of complex poles given by 

dh _|_ 2mm d s ^ Hum ^ ^ 

m d,„ d„, In I 2 g?„. In / 



precisely as in (36). These complex poles have been identified with complex dimensions for fractals [26U27) . motivated 
by the form of the small t asymptotics in (35 1. 

Interestingly, the pole of (d( s ) at s = 1/2 (coming from the (r(2s) factor) has zero residue for all diamonds, and 
so does not contribute to the short time behavior of Kr>(t). Remarkably, this vanishing of the residue at s = l/d w 
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also applies to the analogous zeta function on the Sierpinski gasket [57]. It is believed that this is generic, but this 
has not been proved. The pole of T(s) at s = gives a constant contribution, Cd(0), to Ku{t). The tower of complex 

K(t)/A"i eading (t) 
1.005 




0.001 



FIG. 6: Heat kernel Ko{t) at small time, normalized by the leading non-oscillating term, for the fractal diamond 1)4,2, plotted 
on a log scale in t. The solid [blue] curve is exact; the dashed [red] curve is the approximate expression (431. At such small t, 
these curves are indistinguishable. The relative amplitude of the oscillations remains constant as t — s> 0. 



poles in (42), leads to the oscillatory behavior: 

'lA h -\ 



K D (t) ~ (^^) ^ (ao + 2Re (a^ 2 ^ -'>)) + Cd(0) + - 



where the constants are given explicitly by 



r(s m )C fl (2s m ) 

■7r2s„, 



(43) 



(44) 



This has precisely the form in (34 35 ) , but now in this case of diamond fractals we have explicit expressions for all the 
constants. The leading term cx t~ a/ is multiplied by a periodic function of the form a\ r cos(ln£ sl< ) + an sin(lnt Sli ), 
where a\ r ^ are respectively the real and imaginary parts of <zi, and Sn = 2ir / \xil dw . These oscillatory terms are small 



because of the exponentially small amplitudes of the Gamma functions in ( 44 1 . In Figure [6] we compare this leading 
behaviour with the exact heat kernel trace [dividing out the overall l/t d '^ factor] and find excellent agreement. 
The plot is on a log scale in t, which makes very clear the log periodic oscillations. This figure also shows that the 
magnitude of the log periodic oscillations is very small compared to the leading t dependence. Precisely analogous 
behavior has been found numerically for the Sicrpinksi gasket |49j . 



B. Intuitive physical explanation of complex poles and log periodic oscillations 

A simple intuitive argument for these novel features of the zeta function £(s) and the heat kernel trace K(t) for 
fractals is based on the observation that on fractals the degeneracies of eigenvalues of the Laplacian, and the typical 
gaps in the spectrum of the Laplacian, both tend to grow exponentially 1 rather than polynomially, as they do on a 
regular smooth manifold. This is particularly clear on highly symmetric fractals such as the Sierpinski gasket and the 
diamond fractals, where the eigenfunctions tend to be highly localized and of high degeneracy [33]. Contrast this, for 
example, with the situation on the sphere S d , where the eigenvalues of the Laplacian are 

> j 1 \ j (2n + d- l)(n + d-2)! d _ x , Ar . 
A„ = n(n + d - 1) , deg„ = ^ ( d _ i)i n ! ~ n (45) 

and the zeta function has a meromorphic continuation with simple real poles, the largest of which occurs at s = d/2 
[50] , Correspondingly, the heat kernel trace has the usual Weyl behaviour in 
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FIG. 7: Plot of the ratio of the heat kernel trace function in (51 1, with exponentially growing degeneracies and eigenvalues, 
normalized by its leading small t term. The plots are logarithmic as a function of t, and we have chosen parameters a = 2 
and 6 = 3. The blue plot shows the function (51 1 normalized by the leading term in (52 1, while the red plot includes the first 
log periodic oscillation terms from ( |52[ |. The agreement is excellent at small t. The log periodic oscillations about the leading 
behaviour are clearly visible at small t, and we also see that these oscillations are of very small amplitude, just as for the 
diamond fractals shown in Figure [6] and for the Sierpinski gasket studied numerically in [49] . 



To contrast the consequences of these different large n growths, consider the following "zeta function" with general 
polynomial growth of degeneracies and eigenvalues, n a and n b , respectively, with a and b real parameters. [For the 
Laplacian on a smooth (i-dimensional manifold we would have a ~ d— 1, and b ~ 2, as for the sphere in ( |45[ ).] Then 
the zeta function is 

°° n a 

Cpolynomial(s) = ^ 7^6)7 = Ch(& s ~ <*) , (46) 

n— 1 ^ ' 

which has a simple pole at s — ^ a ~l^ [which reduces to d/2]. The corresponding heat kernel trace function is 

oo 

^polynomial (t) = £ U a e^" 4 . (47) 

ra=l 

A simple Euler-MacLaurin argument yields the small t asymptotics: 

-^polynomial (0 ~ £( a +l)/6 ' 

in agreement with ^ and Efy. 

Now consider the analogous spectral functions with exponentially growing degeneracies and eigenvalues, a n and b n , 
respectively, with a and b real parameters. The zeta function is 

oo a_ 

Ccxponcntial («) = J]] (Un\s = 1 a_ 1 (^9) 

n=l 1 J 6 s 

This zeta function has a tower of complex poles at 

In a 2-Ki . 

s -» = r7 + n > meZ (50) 

In o In 6 

just as for the diamond and Sierpinski gasket fractals. Correspondingly, the heat kernel trace function 

oo 

if exponential (*) = fl " eXp ^ ( 51 ) 
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FIG. 8: Plots of the "Weierstrass function" Schrodinger kernel S(t) defined in (531, with parameters a = 1/2 and b = 3. The 



blue dotted curve is the exact expression, while the black solid line in plots (a), (b), (c) and (d) represent the expression in 
( |54[ | with the oscillatory k sum up to \k\ equal to 0, 1, 10 and 50, respectively. One sees clearly that the agreement is better 
at smaller values of t, and increasing the number of log periodic terms from the k sum improves the agreement out to larger 
values of t. This illustrates the type of wildly oscillatory behaviour expected in a small t expansion of the Schrodinger kernel 
for fractals. 



has log periodic oscillatory behaviour 

l/ln6 



-^exponential \J*) 



fin a/ In b 
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In a 2iri 




2ni , 




r 


lnb_ 


+ 2Re (V 




exp 
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In b In & 


"E6 ln< _ 





(52) 



In FigureP7we compare this small t asymptotics with the function -^exponential (t), showing clear log periodic oscillations 
and excellent agreement with this leading form. As for the diamonds and Sierpinski gasket, the contribution from 
higher complex poles is suppressed by exponentially small Gamma function factors. 

This example makes it clear that the physical origin of the log periodic oscillations can be traced to the fact that 
the degeneracy and eigenvalues grow exponentially with n, rather than as a power of n, as they do on a smooth 
Ricmannian manifold. 



In the case of the Schrodinger kernel ( 18 ) we obtain a highly oscillatory function, whose real part is 



^exponential [tj 



n=0 



a" cos (&" t) 



(53) 



which is precisely the Weierstrass function (we include the n = term, by convention). This function appears in the 
theory of Levy nights and fractals, and the Weierstrass function shows log periodic behaviour [3T]. This function can 
be re- written as 131 



^exponential (0 — ^ ^ 



(-t 2 y 



n=0 K ' v ; 
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which clearly shows the log periodic oscillations. This is difficult to plot when a > 1 because the oscillations are so fast, 
but we can see in Figure [8] the situation for a = 1/2 and b — 3. Keeping an increasing number of log periodic terms 
[corresponding to an increasing number of complex poles of the zeta function], we see better and better agreement 
with the exact function out to larger and larger values of t. The fractal features are self-evident from the plots. 

To conclude this section, we note that this discussion has focused on highly symmetric fractals, those with exact 
self-similarity, motivated by the fact that these fractals exhibit novel features such as log periodic oscillations in the 
heat kernel trace, and associated complex poles of the zeta function. These novel features make them particularly 
interesting from the physical perspective. But it is important to appreciate that there are other fractals, for example 
when randomness or disorder is introduced, where the oscillatory behavior does not necessarily occur. Some specific 
examples, both mathematical and physical, appear in |45[ I51j. 



V. THERMODYNAMICS ON FRACTALS 



Motivated by the problem of defining quantum field theory and quantum gravity on a fractal, we consider first 
a simpler problem: can one define quantum statistical mechanics and thermodynamics on a fractal space? More 
precisely, what is the analogue on a fractal of the usual thermodynamic equation of state 

PV=~ , (55) 

for a smooth d-dimensional closed manifold of volume V, with pressure P and internal energy U? It is not immediately 
clear on a fractal space which dimension plays the role of d, or what is the appropriate volume V, because naively the 
volume of a fractal is infinite. As mentioned already in the introduction, Mark Kac notes in his famous paper "Can 
One Hear the Shape of A Drum?" [34] , that an earlier version of this question had been posed by Lorentz for a regular 
manifold, who asked why the Jeans law of thermodynamics only refers to the total volume of the region occupied by 
a gas, and not to its shape. Here we pose the analogous question for a fractal spatial region. The physical idea is to 
use propagating particles, as they achieve thermal equilibrium, to probe the geometry of the region to which they are 
confined. 

Appropriately for this article, the most efficient way to approach this problem was pioneered by Dowker and 
Kennedy [3J, who studied quantum field theory on spacetimes of the form S 1 x M., with S 1 being the compactified 
Euclidean time of period fij3 — Ti/T for the Matusbara modes. To illustrate this idea consider a massless bosonic 
field [a "photon"] confined to live on a spatial region that is a fractal T . The thermodynamic partition function is 
expressed as 

lnZ(T, V) = ~ lnDet 51x ^ (d 2 /dt 2 - c 2 A) (56) 
where A is the Laplacian on the spatial region J 7 , and c is a velocity parameter [the "speed of light"]. The expression 



( 56 1 is usually derived using a mode decomposition of the bosonic field, summing over states in phase space, which 
in turn relies on an analysis of Fourier transforms. But on a fractal there is, as yet, no well-defined notion of Fourier 
analysis. However, we can by-pass this step as follows. Recall that the spectral partition function for a single bosonic 
oscillator of frequency u) is 

lnZ(J» = _^-ln(l-e-^) 



3G 



2 ^ n 

71=1 



(3Tiuj ^ f3h f°° dr 



E 



2 ^ y^TT Jo T3/ 2 



e e 



1 r°° dr 

-e 



2 Jo r 



r e~W T (57) 



n— — oo 



where in the last step we have used the Poisson summation formula. Using the Kubo, Martin and Schwinger (KMS) 
condition [S3J for thermodynamic equilibrium of a bosonic field at temperature T, we recognize the second factor 
in the integrand as a sum over Matsubara frequencies (^f ) 2 , corresponding to the spectrum of the operator d 2 /dt 2 
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with periodic boundary conditions <f>(t + h/3) = <j)(t). Identifying ur = c k with the eigenvalues of — c A, and tracing 



over all modes [62) , we recover expression ( 56 1 , which can be written as follows 



\nZ(T,V) 



dr 



6 3 0,e 



KAr) 



(58) 



where O3 is the Jacobi theta function and we have defined the "photon" thermal wavelength, Lp = j3hc. All 
information about the spatial region J-, including possible dependence on its dimension and volume, is contained in 
the heat kernel trace factor Kj?(t), while the thermodynamic information is encoded in the theta function factor. 
The form in (58) is particularly convenient because it factors, in the integrand, the temperature dependence from the 



spectral dependence of the Laplacian on the spatial region. 

The heat kernel Kjt(t) contains an implicit length scale L s , which we call the "spectral length", that defines a 
dimensionless Laplacian A = L 2 S A [recall that the eigenvalues of A have dimensions of 1/length 2 ]. Results in the 
mathematical literature refer to the dimensionless Laplacian A. For black-body radiation confined to manifolds of 
complicated shape, it is difficult to make an explicit mode decomposition and find an explicit expression for the heat 
kernel trace, except for highly symmetric regions such as hyper-cubes or spheres and sections of spheres [54] . However, 
we can learn about the thermodynamic [large volume] limit from the Weyl expansion of the heat kernel trace. Note 
that the large volume thermodynamic limit correspon ds to L s 3> Lp, which is a high temperature limit T 3> hc/L s . 
This probes the small r behavior of K.jt{t / L 2 S ). Using (35), the leading term is 



cLI 



( 47rr )d/2 



(59) 



which we can use to define the "spectral volume" V s [55] 

V s = cL d J 



where the numerical coefficient is determined below [see (66 



(60) 

Physically, we are using the photons in thermal 
1 they are confined. This physical idea gives a direct 



equilibrium to define the volume of the spatial region in whic 
link between the thermodynamic limit and the leading term of the Weyl expansion: 

thermodynamic limit O leading term of Weyl expansion 



(61) 



For photons confined to a regular manifold region, this reproduces all the standard results of thermodynamics, 
including the equation of state (55). To see this, recall the thermodynamic limit following from the leading Weyl 
behaviour Q: 



In Z 
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This leads immediately to the usual equation of state ( 55 ) 



We can also formulate this in terms of the zeta function £(s) (see again the work of Dowker and Kennedy [3]) 
a regular manifold Ai, straightforward manipulations show that the partition function can be expressed as 



lnZ(T, V) = - 



L 



C 



M 



1 

7TC 



L 



2 s 



r(2a) 0?(2s + 1) Cm(s) ds 



(62) 



For 



(63) 



The first term gives the standard zero temperature "vacuum energy" contribution [3] , proportional to ^4 (— |), 
which has recently been generalized to quantum graphs [56] . The second term in ( |63[ ) encodes finite temperature 
corrections to the internal energy U; the Riemann zeta factor Cr(2s + 1) arises from the sum over Matsubara modes. 
The large volume behaviour is controlled by the largest pole of Cm( s )i which occurs at s — d/2, and reproduces the 
thermodynamic limit in ( 62 ) , and also determines the volume V in terms of the zeta function: 



V = (47r) d ^T(d/2)L d Re S [C M (s)] s= 



d/2 



(64) 
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On a fractal spatial region T the same procedure follows through, replacing (m( s ) with Cf( s )- Together with the 
definition of the spectral volume in ( 60 1 , the equation of state is found to be |55| 



PV S = -*j-U 



(65) 



This has the same form as (55), but physically is quite different. First, note that the spectral volume scales with the 
length L s to the power d s , the spectral dimension, not the Hausdorff dimension. Consequently, the dimension d in 



(55) is replaced by the spectral dimension d s in the fractal equation of state (65). The numerical coefficient in the 



spectral volume (60) is fixed by the residue of the zeta function at the pole at s = d s /2: 



V s = (4^/ 2 r(d s /2)LfRes[^( S )] s=( ^ 



/2 



(66) 



Given that there is no well-defined notion of Fourier transforms on fractals, it is somewhat surprising that this 
thermodynamic analysis is possible. A standard approach to thermodynamic equilibrium is in terms of momentum 
modes, as we write 



InZ 



V m 



d d k 
{2ir) d 



flUJ 

~2~ 



1 



In (1 - 



(67) 



where we have explicitly written V mo m for the momentum space volume. On a regular smooth manifold V mom ~ 1/V, 
the inverse of the spatial volume, and so V enters the thermodynamic equation of state. On the other hand, it is 
clear that the V in the equation of state is really 1/Vmom, and on a fractal this is replaced by the spectral volume V s 
defined in (66). This suggests a provocative idea, which deserves to be explored further: on a fractal-like region, it 



appears that the spectral volume (which we might identify with momentum space) has, in some appropriate sense, 
a dimension of d s , while the spatial volume has a dimension of d^. On a smooth manifold these dimensions are 
the same, but on a fractal they are generically different. This could lead to some potentially unusual properties of 
thermodynamic equilibrium, which deserve to be probed experimentally. 



VI. APPROXIMATE FRACTALS 



The possibility of probing experimentally such strange properties of fractal regions raises an immediate objection 
that one can never fabricate an infinitely iterated fractal such as a perfect Sierpinski gasket or diamond, or even a 
Sierpinski carpet [see Figure [9] , because there is always a lower limit on the spatial scale that can be realized in a 
real physical system. However, it turns out that these strange properties [for example, the spectral dimension and 
log periodic oscillations] can be seen already with relatively modest level of iteration of the fractal. It is a question of 
resolution scale, and can be illustrated as follows. 



FIG. 9: A Sierpinksi carpet, 
numerically, estimated as d 3 ; 



Its Hausdorff dimension is dh = In 8/ In 3 « 1.89, but its spectral dimension d s is only known 
i 1.80525, and rigorously bounded as 1.674 < 21n8/lnl2 < d s < 2 In 8/ ln(28/3) < 1.86 [57]. 



To be very explicit, we consider again the diamond fractals. For a finite level of iteration, a diamond fractal is 
undoubtedly a one-dimensional quantum graph. Therefore, in the asymptotic small t limit, the heat kernel trace 
should behave in a way governed by the one-dimensional Weyl law (pi, rather than by the fractal Weyl law (35 ). This 



is shown in Figure 10 where we plot the ratio of the heat kernel trace K(t) to the leading fractal value from (34), for 
5 and 10 levels of iteration of the diamond. This is plotted as a function of log t, for several different ranges, probing 
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FIG. 10: Plots of the heat kernel trace K{t), normalized by the leading fractal form in ( 35 1, for the Dg t 2 diamond fractal, plotted 
on a log scale in t. The solid [blue] curves show the result for 10 iterations of the fractal, while the dotted [red] curves show the 
result for 5 iterations. The dot-dashed [black] curves show the expected result for a one-dimensional quantum graph. The first 
plot probes t down to 10 -6 , while the second probes t down to 10 -8 . Notice the region of t in which these approximate fractals 
behave as if they had spectral dimension d s — In 6/ In 2 « 2.58, and also exhibit log periodic oscillations with the expected 
period. Also note that this region extends down to smaller values of t as the number of iterations increases. 



smaller and smaller values of t. We also show with the dashed line, the one-dimensional behaviour obtained from the 
quantum graph one-dimensional limit in (JsJ) . These plots show clearly that at very small t, with a finite number of 
iterations of the fractal, the heat kernel trace indicates truly one-dimensional behaviour. As the number of iterations 
increases, we see a decrease in the value of t at which this one-dimensional behaviour sets in. This is to be expected. 
But the surprising thing is that there is a well-defined range of t in which the finitely-iterated "approximate fractal" 
behaves as if it has spectral dimension d s . Not only does d s /2 [rather than 1/2] give the correct overall scaling of the 
heat kernel trace, we also see the log periodic oscillations, indicative of "would-be" complex poles in the zeta function, 
even though the approximate fractal is a one-dimensional graph. The scale at which this occurs is determined by 
the size of the smallest length scale of the fractal at that level of iteration. Furthermore, this is clearly visible, for 
some window of t, with only 5 fractal iterations. This suggests that a physical probe of these fractal properties might 
be possible by choosing a propagating probe of the appropriate wavelength to match this scale, in a region in which 
the approximate fractal looks like an ideal fractal. Similar results have been obtained [JO, 58 , both numerically and 
experimentally, for Fibonacci stacks of dielectrics [5], which are spatially regular but spectrally fractal. This general 
question of when an approximate fractal behaves like an ideal fractal is a crucial physical question, and it deserves 
further mathematical attention, both analytic and numerical, in order to be made more precise. 
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VII. OFF-DIAGONAL HEAT KERNEL PROPAGATOR 



Another vitally important problem for which there are only partial mathematical results for fractals is the question 
of how the well-known heat kernel propagator in (13 14) generalizes to fractals. Here the main rigorous result [2lTf2"3] 



gives upper and lower bounds indicating sub-diffusive propagation: 



Cl 
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(68) 



for some real positive constants ci , ...,04, with d w > 2 being the walk dimension and p(x, y) an appropriate geodesic 
distance between x and y. Very little is known about how P t (x,y) behaves between these bounds, although it can 

be shown to be a log periodic function of the ratio p(x, y) d ^~ 1 jt d ^- 1 [591 EQ]- Taking the functional form of one 
of these bounds as an averaged asymptotic form for P t (x,y), we learn something interesting about the form of the 
Green's function at short distances. Recalling the relation (15) between G(x,y) and P t (x,y) = K(x,y;t), we deduce 
the asymptotic form 



G(x,y) 



( P (z,y))M<*s-2)/2 (p( X ,y))dH 



d h -d u 



d s ^2 



(69) 



Thus, the familiar Green's function exponent (d — 2) in ( |16[ ) for a smooth d-dimensional manifold, is replaced on a 
fractal by (dh — d w ), or equivalently by d w (d s — 2)/2. This shows that for propagation of diffusion on a fractal the 
critical dimension is when the spectral dimension d s [not the Hausdorff dimension dh] is equal to 2. Implications for 
statistical mechanics and Bose-Einstein condensation have been studied recently by Chen [51]. In this critical case, 
we would expect logarithmic behaviour instead of (69). 



VIII. CONCLUSION 



In conclusion, the zeta function and heat kernel continue to be useful spectral functions in describing physical 
particles or fields propagating on fractal spaces. However, they exhibit some interesting new features on fractals, 
quite different from their behaviour on smooth manifolds, and these features lead to new mathematics and potentially 
interesting new physics. 
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